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ABSTRACT 

Our recently completed one billion particle Via Lactea II simulation of a Milky Way-sized dark 
matter halo resolves over 50,000 gravitationally bound clumps orbiting today within the virialized 
region of the main host. About 2,300 of these subhalos have one or more "progenitors" with M > 
10 6 M Q at redshift z = 11, i.e. massive enough for their gas to have cooled via excitation of H2 and 
fragmented prior to the epoch of cosmic reionization. We count 4,500 such progenitors: if these were 
able to convert a fraction of their gas content into very metal-poor stars with a Salpetcr initial mass 
function (IMF), they would be shining today with a visual magnitude My = 6.7 per solar mass in stars. 
Assuming a universal baryon fraction, we show that mean star formation efficiencies as low as 0.1% 
in progenitors <C 10 s M© would overproduce the abundance of the faint Galatic dwarf spheroidals 
observed by the Sloan Digital Sky Survey. Star formation at first light must have occurred either with 
an IMF lacking stars below 0.9 M©, or was intrinsically very inefficient in small dark matter halos. If 
the latter, our results may be viewed as another hint of a minimum scale in galaxy formation. 
Subject headings: cosmology: theory - dark matter - galaxies: dwarfs - halos 



1. INTRODUCTION 

Galaxy halos are one of the crucial testing grounds for 
structure formation scenarios, as they contain the im- 
prints of past accretion events, from before the epoch of 
reionization to the present. In the standard ACDM con- 
cordance cosmology, objects like the halo of our Milky 
Way are assembled via the hierarchical merging and ac- 
cretion of many smaller progenitors. Subunits collapse at 
high redshift, are dense, and have cuspy density profiles, 
and when they merge into larger hosts they are able to 
resist tidal disruption. Indeed, the Via Lactea Project, 
a suite of some of the largest cosmological simulations 
of Galactic dark matter substructure, has shown that 
galaxy halos today are teeming with surviving "subha- 
los" (Diemand et al. 2007, 2008). The predicted sub- 
halo counts vastly exceed the number of known satel- 
lites of the Milky Way, a "substructure problem" that 
has been the subject of many recent studies. While a 
full characterization of this discrepancy is hampered by 
luminosity bias in the observed satellite luminosity func- 
tion (LF) (Koposov et al. 2008), it is generally agreed 
that cosmic reionization may offer a plausible solution 
to the apparent conflict between the Galaxy's relatively 
smooth stellar halo and the extremely clumpy cold dark 
matter distribution. In this hypothesis, photoionization 
heating after reionization breakthrough reduces the star 
forming ability of newly forming halos that are not suffi- 
ciently massive to accrete intergalactic gas (e.g. Bullock 
et al. 2000; Somerville 2002; Kravtsov et al. 2004; Ri- 
cotti & Gnedin 2005; Moore et al. 2006; Strigari et al. 
2007; Simon & Geha 2007; Madau et al. 2008). The ex- 

1 Department of Astronomy & Astrophysics, University of 
California, Santa Cruz, CA 95064, pmadau, dicmand@ucolick.org. 

2 Institute for Advanced Study, Einstein Drive, Princeton, NJ 
08540, mqk@ias.edu. 

3 Institute for Theoretical Physics, University Zurich, 
Winterthurerstrasse 190, 8057 Zurich, Switzerland, dpot- 
ter,moore,stadel@physik.uzh.ch. 

4 Astronomy Department, University of Michigan, Ann Arbor, 
MI, 48109, mzemp@umich.edu. 



act value of this mass threshold remains uncertain, as it 
depends on poorly known quantities like the amplitude 
and spectrum of the ionizing background and the abil- 
ity of the system to self-shield against UV radiation (e.g. 
Efstathiou 1992; Dijkstra et al. 2004). 

In this Letter we take a different approach and focus 
instead on the sources of reionization itself. We use our 
one billion particle Via Lactea II simulation to constrain 
the character of star formation (initial mass function and 
efficiency of gas conversion into stars) at first light, i.e. in 
subgalactic halos prior to the epoch of reionization. The 
starting point of our investigation is the finding that as 
many as 2,300 bound clumps that survive today in the 
Via Lactea II halo have early progenitors that were mas- 
sive enough for their gas to cool via excitation of molec- 
ular hydrogen and fragment before reionization break- 
through. If low-mass stars formed in the process, their 
hosts would be shining today as faint Milky Way satel- 
lites. 

2. VIA LACTEA II 

Our recently completed Via Lactea II simulation, one 
of the highest-precision calculation of the assembly of the 
Galactic CDM halo to date (Diemand et al. 2008), offers 
the best opportunity for a systematic investigation of the 
fossil records of reionization in the halo of the Milky Way. 
Via Lactea //employs just over one billion 4, 100 M© par- 
ticles to model the formation of a M200 = 1-9 x 10 12 M© 
Milky- Way size halo and its substructure. It resolves 
50,000 subhalos today within the host's r2oo = 402 kpc 
(the radius enclosing an average density 200 times the 
mean matter value). 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
5-year data require the universe to be fully reionized by 
redshift z = 11.0 ± 1.4 (Dunkley et al. 2008). Assuming 
that the region around the Milky Way was reionized at 
about the same epoch of the universe as a whole, we can 
trace the progenitors of present-day surviving substruc- 
ture back to a time prior to reionization breakthrough, i.e. 
before star formation was quenched by photohcating. We 
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dwarf galaxy of a descendant subhalo. We include in our 
analysis only progenitors that either survive individually 
during the hierarchical clustering process or contribute 
to a gravitationally-bound descendant at z = 0. About 
10,000 "first generation" systems above 10 6 M© are to- 
tally disrupted by tidal forces: for half of them all of 
their particles lie today within r2oo and contribute to the 
smooth stellar halo. Figure [1] shows an image of all Via 
Lactea //progenitor halos above 10 6 M© at z = 11 with 
a bound descendants within ?'2oo today, while Figure [2] 
shows the cumulative mass functions of progenitors and 
parents. We track forward in time 4,500 of these early 
minihalos into 2,300 descendants today. 5 
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Fig. 1. — Top: Projected dark matter density map of Via 
Lactea II at z = 11. The image covers 3.5 comoving Mpc across 
and shows (magenta) all progenitor minihalos with M > 10 s Mq 
that have a bound descendant today within r200- Bottom: the 
2 = descendants of those early progenitors. The image covers 
800 kpc across. 

use the 6DFOF group finder described in Dicmand ct al. 
(2007) to identify peaks in phase-space density at z = 11 
and z = 0. The resulting groups contain between 16 
and a few thousand particles linked together in the cen- 
ters of halos and subhalos. Note that the gravitationally 
bound region often extends beyond this central group. 
We link a descendant z = group "A" to a high rcdshift 
progenitor "B" if "A" contains more than 10 particles 
from "B" and more than any other descendant of "B". 
Thus a descendant can have more than one progenitor, 
but a progenitor is linked to at most one descendant. 
Using the halo centers provided by 6DFOF ensures that 
material from "B" does indeed contribute to the central 
regions of "A" , i.e. that a star cluster that formed in the 
core of a progenitor clump would end up today in the 



Fig. 2. — Via Lactea II cumulative substructure mass function 
at different epochs. Solid curve: All 4,500 progenitor halos with 
M > 10 6 Mq at z = 11 that contribute to a gravitationally bound 
descendant (2,300 of them) today within T2Q0- Dot-dashed curve: 
same for the 2,900 M > 10 6 Mq progenitors at z = 13. Also shown 
for comparison is the mass function of all Via Lactea II subhalos 
at the present epoch {dashed curve). 

3. STAR FORMATION AT FIRST LIGHT 

Cosmological hydrodynamics simulations of structure 
formation in a ACDM universe have found that in the 
early collapse of > 10 6 M Q systems, enough H 2 is pro- 
duced to cool the gas and allow star formation within a 
Hubble time (e.g. Abel et al. 2002; Bromm et al. 2002). 
Figure [5] shows that a few hundred progenitors are above 
the "atomic cooling" mass threshold of - 3 x 10 7 M© 
(Tvir ^ 10 4 K), where gas can cool and fragment via ex- 
citation of hydrogen Lya. Numerical studies also suggest 
that, while primordial stars were likely very massive and 
formed in isolation at z > 20 (e.g. O'Shea & Norman 
2007), low-mass second- generation stars can form as soon 
as a minimum pre-enrichment level of Z = 10~ 5±1 Zq is 
reached (e.g. Schneider et al. 2002). The existence of 

5 Note that the size of the Galaxy halo is not known very pre- 
cisely, since the effects of disk formation on the total mass distri- 
bution are still poorly constrained. When comparing Via Lactea II 
to the Milky Way, this translates into an uncertainty in the abun- 
dance of subhalos of about a factor of two (Klypin et al. 2002; 
Dutton et al. 2007). 
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0.8 M Q extremely iron-deficient stars in the halo of the 
Milky Way (Christlieb et al. 2002; Frebel et al. 2005) 
indicates that low-mass star formation was possible at 
very- low mctallicities (e.g. Tumlinson 2007). In or- 
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Fig. 3. — The luminous fossil remnants of reionization. Solid 
curve: present-day LF of all Via Lactea II subhalos that had 
one or more progenitor with M > 10 6 Mg at z = 11. At this 
epoch, a fraction /* of the gas content of each progenitor is turned 
instantaneously into stars with a Salpeter IMF and metallicity 
Z = Z /2OO. Dashed curve: same for M > 5 X 10 5 M . Dot- 
dashed curve: same for a Chabricr (2003) IMF. Dotted curve: same 
for Z = Zq/5. In all the above cases visual magnitudes at age 13.7 
Gyr have been calculated using Bruzual & Chariot's (2003) models 
(1994 Padova tracks). 

der to constrain the number of low-mass stars that lit 
up dark matter halos before reionization, our default 
model makes a number of simplifying assumptions: 1) 
all progenitor hosts above M(z = 11) = 10 6 M© have 
a gas content given by the universal baryon fraction, 
M gas = (tt b /tt m )M = 0.17 M; 2) at z > 11 a fraction 
/» of this gas is turned into very metal-poor stars with 
Z = Z Q /200 (the lowest metallicity allowed by Bruzual 
& Chariot's 2003 stellar population synthesis models) 
and a standard Salpeter initial mass function (IMF) in 
the range 0.1 < m* < 100 M . The star formation ef- 
ficiency is independent of M; 3) star formation is sup- 
pressed at later epochs in all progenitors and their de- 
scendants; 4) primordial stellar systems are deeply em- 
bedded in progenitor minihalos and remain largely un- 
affected by tidal stripping even if their hosts are not. 
The complete tidal disruption of a host, however, also 
destroys its stellar system. 

We are now in the position to construct the present-day 
LF of the fossil remnants of reionization in the halo of 
the Milky Way. According to Bruzual & Chariot (2003), 
a stellar population undergoing an instantaneous burst 
will be shining at age 13.7 Gyr with a visual magnitude 
per solar mass in stars of 6.7 mag (Salpeter IMF, Z = 
Zq/200). The steep LF of Via Lactea II first-generation 
halos predicted in this case is shown in Figurc[3]for (/*) = 
0.03. Luminous remnants are distributed between My = 
—2 and My = — 10, and there are more than two hundred 



relatively bright objects with My < —6, i.e. as bright 
as Bootes. A mean star formation efficiency three times 
smaller would shift this curve 1.2 mag to the left. The 
cutoff at faint magnitudes simply reflects our assumption 
that clumps below M{z = 11) = 10 6 M© never form 
stars: the effect of lowering this mass threshold by a 
factor of two is also shown for comparison. Note that 
implausibly increasing the metallicity of early stars to 
Z = Zq/5 would only shift the curve 0.3 mag to the left 
(Bruzual & Chariot 2003). Similarly, a Chabrier (2003) 
IMF in the same mass range instead of Salpeter has only 
a small (0.2 mag) brightening effect. 

4. CONSTRAINTS FROM THE SDSS 

Over the last two years, the Sloan Digital Sky Sur- 
vey (SDSS) has doubled the number of known Milky 
Way dwarf spheroidals (dSphs) brighter than My = — 2 
(e.g. Bclokurov et al. 2007). An accurate estimate of 
the total number of dwarfs requires a correction to the 
observed LF that accounts for completeness limits and 
that depends on the unknown intrinsic spatial distribu- 
tion of sources. Rather than producing an "unbiased" 
list of Galactic satellites as in Tollerud et al. (2008), we 
apply here the completeness limit of the SDSS DR5 to 
our predicted fossil remnants in Via Lactea II (Fig. [3]) 
to construct an artificial DR5 sample of luminous first- 
generation substructure. As computed by Koposov et al. 
(2008) (and fitted by Tollerud et al. 2008), the maximum 
accessible distance beyond which an object of magnitude 
My will go undetected by the SDSS is 
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where /drs = 0.194 is the fraction of the sky covered by 
DR5. Therefore, only subhalos brighter than My — —6.6 
will be detected all the way to t , 2oo : first-generation rem- 
nants with My = —4 will be included in the sample only 
if they are closer than 120 kpc. (The above detection 
limit applies to objects with surface brightness < 30 mag 
arcsec -2 , i.e. we are implicitly assuming that our lumi- 
nous progenitors are less diffuse than this.) As a first or- 
der correction we also assume that subhalos within r max 
are detected with 100% efficiency, and account for the 
partial sky coverage of SDSS by multiplying the number 
of detectable subhalos by /drs ■ The resulting cumulative 
LF N DR5 (< My) is plotted in FigureHfor values of (/*) 
that decrease from 3% down to 3 x 10~ 4 . In this range, 
the faint-end of the predicted LF is dominated by the 
remnants of 10 6 — 10 7 M© progenitors. For an efficiency 
of 1%, SDSS should have detected ~ 150 first-generation 
systems brighter than My = —2, while only a dozen 
are actually observed. The known Milky Way satellites 
are also plotted for comparison: these include the SDSS 
eleven dwarfs with My < — 2 (including Leo T which, at 
a distance of 417 kpc, lies just outside r2oo) as well as the 
eleven (times /drs to correct for the partial sky cover- 
age) classical (prc-SDSS) bright satellites. Clearly, mean 
star formation efficiencies as low as 0.1% in progenitors 
below 10 7 Mq would overproduce the abundance of ultra- 
faint dSphs. The mean efficiency of star formation must 
actually scale with halo mass in order to reproduce the 
observations with first-generation remnants. This is be- 
cause a constant efficiency translates into a steep LF that 
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follows the steep mass function of CDM field halos. By 
contrast, the 23 known satellites of the Milky Way have 
a flat LF and shine with luminosities ranging from about 
a thousand to a billion times solar. The dashed curve 
in Figure |4] shows a toy model with a mass-dependent 
efficiency (/*) = (0.02,0.0025,0) for M in the range 
M = (> 7 x 10 7 ,3.5 x 10 7 -7 x 10 7 , < 3.5 x 10 7 M ). 
Note how the data allow efficiencies of order 1% only 
above the atomic cooling threshold. 
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Fig. 4. — Via Lactea II first-generation (z = 11) subhalos de- 
tectable today by the SDSS DR5. All curves assume M > 10 6 M Q , 
a Salpeter IMF, and Z = Zq/200. The data points are taken 
from the compilation of Tollerud et al. (2008). Labelled solid 
curves: (/») = 0.03,0.01,0.003,0.001,0.0003. Dashed curve: mass- 
dependent star formation efficiency, (/») = (0.02,0.0025,0) for M 
in the range (> 7 X 10 7 ,3.5 X 10 7 -7 X 10 7 , < 3.5 X 10 7 M Q ). Also 
plotted for comparison is the LF of the remnants of z = 13 halos 
with (/*) = 0.0005 (short dot-dashed curve at the bottom). 



5. SUMMARY 

In this work we have used results from the one billion 
particle Via Lactea II simulation to study for the first 
time the fossil signatures of the pre-reionization epoch 
in the Galactic halo and set constraints on the baryonic 
building blocks of today's galaxies. We have traced the 
progenitors of present-day surviving substructure back 
to a time prior to reionization breakthrough, i.e. before 
star formation was quenched by an external UV back- 
ground. We have then populated early-forming mini- 
halos with very metal-poor stars following a Salpeter 
IMF, and looked at the photometric properties today 
of their 2,300 descendants within 7"2oo. We have shown 
that even star formation efficiencies as low as 0.1% in 
progenitor minihalos below 10 M would overproduce 
the abundance of Galatic dSphs observed by the SDSS. 
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This value should be regarded as an upper limit since 
luminous dSphs are known to have formed a fraction of 
their stars after reionization (e.g. Orban et al. 2008). It 
is also meant to be taken as a mean value, since an alter- 
native possibility is to allow stars to form efficiently in 
a small fraction of minihalos at some given mass and to 
suppress star formation in the remainders, rather than 
reducing the efficiency across the board. Note that the 
parameter /* is degenerate with the baryon fraction: if 
the gas content of early virializcd structure was smaller 
than the universal value (as found, e.g., by O'Shea & 
Norman 2007), then our limits on the star formation ef- 
ficiency would be correspondingly higher. 

We can now test a posteriori the main assumption un- 
derlying our work, i.e. that early Via Lactea II pro- 
genitors were populated with "pre-reionization" stellar 
systems. In the simulated high-resolution volume, there 
are 17,700 M > 10 6 M subhalos at z = 11, spread 
over a comoving volume of 78 Mpc 3 containing a total 
mass of Mdm = 3.5 x 10 12 M Q . Stars distributed ac- 
cording to a Salpeter IMF produce during their lifetime 
/ 7 w 4, 000 Lyman continuum photons per stellar proton, 
of which only a fraction / esc will escape into the inter- 
galactic medium (IGM). The total stellar mass in these 
early-forming progenitors is M* = 3.6 x 10 10 (/») M Q . 
Hydrogen photoionization requires (1 + iV rec ) photons 
above 13.6 eV per atom, where N ICC is the number of 
radiative recombinations over a Hubble time. The to- 
tal mass of intergalactic gas that can be kept ionized is 
then M ion = / 7 M*/ esc /(l + N Iec ). The condition M ion < 
(n b /Q m ) M DM then implies (/») < 0.004 [(l+N Tec )/ f esc ], 
where the factor in square brackets is of order 10 or so. 
The low star formation efficiencies derived in this work 
appear then to be consistent with the idea that the region 
around the Milky Way was reionized at z < 11 by exter- 
nal radiation, before stars in the Local Group formed in 
sufficient numbers (Weinmann et al. 2007). 

We conclude that star formation at first light must 
have occurred either with an IMF lacking stars below 
0.9 Mq, or was intrinsically very inefficient in small dark 
matter halos. If the former, this may be an indication 
of an upward shift in the mass scale of the IMF at early 
cosmological times owing, e.g., to the hotter cosmic mi- 
crowave background (Larson 2005). If the latter, our 
results may be viewed as another hint (see Strigari et 
al. 2008) of a minimum scale in galaxy formation, below 
which supernova feedback (Dekel & Silk 1986) and/or 
H2 photodissociation by a Lyman- Werner background 
(Haiman et al. 1997) sharply suppress star formation. 
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